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ABSTRACT 



The intermediate-frequency peaked BL Lacertae (IBL) object 3C 66A is detected 
during 2007 - 2008 in VHE (very high energy: E > 100 GeV) 7-rays with the VERITAS 
stereoscopic array of imaging atmospheric Cherenkov telescopes. An excess of 1791 
events is detected, corresponding to a significance of 21.2 standard deviations (a), in 
these observations (32.8 hours live time). The observed integral flux above 200 GeV is 
6% of the Crab Nebula's flux and shows evidence for variability on the time-scale of 
days. The measured energy spectrum is characterized by a soft power law with photon 
index T = 4.1 ±0.4 s t a t i0.6 S y S . The radio galaxy 3C 66B is excluded as a possible source 
of the VHE emission. 

Subject headings: galaxies: active — BL Lacertae objects: individual (3C 66A) - 
gamma rays: observations 



1. Introduction 



Wills &: Wills! ()1974l ) first identified 3C 66A as a QSO using optical observations. It was 
subsequently classified as a BL Lac object based on its significant optical and X-ray variability 
(jMaccagni et al.lll987l ). BL Lac objects are characterized by a double- humped spectral energy 
distribution (SED) and are further classified according to the location o f the lower e nergy hump, 



usually interpreted as synchrotron emission from relativistic electrons. iPerri et al.l (|2003l ) locate 
the synchrotron peak between 10 15 and 10 16 Hz. Therefore 3C 66A is classified as an intermediate- 
frequency peaked BL Lac (IBL). To date, the majority of BL Lacs detected at VHE (very high 
energy: E > 100 GeV) are HBLs (hi gh-frequency peake d BL Lacs). Only one other IBL, W Comae, 
has been detected above 100 GeV (jAcciari et al.l 12008). 



During states of high flux, continuum emission from the jet is dominant and overshadows the 
few emission lines from the rest of the galaxy. Blazars have few, if any, detectable emission lines, 
which makes determining the redshift difficult even under the best conditions. Base d on a single 
line, interpreted as Mg I I, 3C 66A was determ ined to be at a redshift of z = 0.444 (jMiller et al 



19781 ). In addition to this. lLanzetta et al.l (|1993l ) identified a weak Ly-alpha line corroborating these 
results. Si nce both measurem ents rely on a single line, the redshift of this BL Lac is considered 
uncertain. iFinke et al.l (|2008l ) recently derived a lower limit of z = 0.096. The determination of the 
redshift is crucial to under standing this source at VHE energies due to the effect of the extragalactic 



background light (EBL) (jHauser &; Dwekl l200ll ). VHE 7 rays are abso rbed via pair productio n 
interactions with the infrared component of the EBL (^ v helebi —* e + e~) (jGould &: Schrederlll967l ). 
At VHE energies, this absorption causes a decrease in the observed flux and a softening of the 
observed spectrum. One can calculate an optical depth (t(z,E)) based on an EBL density model, 
the redshift and the 7-ray energy. The optical depth can be used to calculate the flux corrected 
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for extragalactic absorption from the observed flux at a given energy {Fi n t = e T ^ z ' E ^F f JS ). Without 
an accurate measure of the redshift, an accurate photon spectrum intrinsic to the blazar cannot be 
calculated or modeled. 

The EGRET source 3EG J0222+4253, detected at an integral flux b etween (12.1 ±3.9) x 
1CT 8 and (25.3 ± 5.8) x 1(T 8 photons cm~ 2 s -1 , is associated with 3C66A (jHartman et al.lll999n . 
Measurements of the spectrum indicated that the spectral index of 2.01 ± 0.14 was influenced by 
the nearby pulsar PSR 0218+42, which is also inside the EGRET error box. A detailed study of the 
energy-dependent position of the EGRET source shows that the highest energy photo ns are coming 



from the BL Lac and thus the spectrum is t hought to continue out to the VHE band (IKuiper et al 



2000). A re-analysis of the EGRET data bv lNandikotkur et al.l ([20071 ) finds a harder spectral index 



of 1.95 ± 0.14 and a flux above 100 MeV of (17.7 ± 2.8) x 10 -8 photons cm" 2 s _1 . Recently the 
Fermi Gamma-ray Space Telescope also reported a detection of 3C 66A, contemporaneous with 
VHE data t aken by VE RITAS in the 2008-2009 season, at a higher flux than previously reported 
by EGRET ^Tostil boosl ) . 



There have been several attempts to detect 3C 66A in VHE 7 rays. The Crimean Astrophysical 
Observator y reported a 5.1 a det ection above 900 GeV at an average integral flux o f 2.4 x 10~ ] 1 



cm 2 s 1 (jStepanyan et al.ll2002l ). Additionally, both the Whipple 10 m telescope (jHoran et al 



2004 ) and HEGRA telescope array (jAharonian et al.1 12000 ) observed 3C 66A and reported upper 



1 1 —2 - 

cm s 



limits of < 0.35 x 10~ 11 cm~ 2 s~ 1 (99.9% confidence; above 350 GeV) and < 1.4 x 10" 
confidence; above 630 GeV), respectively. STACEE reported several upper limi ts between < 1.0 an d 



< 1.8x 10 10 cm 2 s 1 above 150 - 200 GeV depending on the source spectrum ( Bramel et alJl2005l ). 



Recently, MAGIC reported a 5.4 a detection of VHE emission above 150 GeV from observations 
in September to December 2007 coincident with 3C66B. T hey exclude 3C 66A as the source of the 
VHE emission at an 85% confidence level (jAliu et al.ll2009l ). 



2. VERITAS Detector & Observations 



The VERITAS detector is an array of four 12 m diam eter imaging atmospheric Cherenkov 
telescopes located in southern Arizona (jWeekes et al.l |2002j) . Designed to detect emission from 
astrophysical objects in the energy range from 100 GeV to greater than 30 TeV, VERITAS has an 
energy resolution of ~15% and an angular resolution (68% containment) of ~ 0.1° per event. A 
source with a flux of 1% of the Crab Nebula is detected in ~ 50 hours of observations while a 5% 
Crab Nebula flux source is detected in ~ 2.5 hours. The field of view of the VERITAS teles copes 
is 3.5°. For more details on the VERITAS instrument and technique, see lHolder et al.l (|2008l ). 



VERITAS observed 3C 66A for 14 hours from September 2007 through January 2008 (hereafter, 
the 2007-2008 season). From September through November 2008 (hereafter, the 2008-2009 season), 
a further 46 hours of data were taken. In total, 180 twenty-minute exposures were made, where 
109 exposures passed selection criteria which remove data with poor weather (based on infrared 
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sky-temperature measurements) or with hardware-related problems. Data collection frequently 
occurred during poor weather conditions causing the lower selection throughput demonstrated 
here. In total, the 2007-2008 season resulted in 4.7 hours live time and the recent 2008-2009 season 
produced 28.1 hours live time. The average zenith angle was 17.3°. All data were taken on moonless 
nights in "wobble" mode where the teles copes are pointed away from the source by ±0.5° to allow 
for simultaneous background estimation (jBerge et al.l 120071 ). 



3. Analysis Methods 



Prior to event selec tion an d bac kgroun d subtraction , the s hower images are calibrated and 
cleaned as described in ICoganl (|2006l ) and iDaniel et al.l (|2007l ). Several noise-reducing event- 
selection cuts are made at this point, including rejecting those events where only the two closest- 
spaced telescopes participated in the trigger. Follo wing the cal ibration and cleaning of the data, the 
events are parametrized using a moment analysis rtHillaslll985T). From this moment analysis, scaled 



param eters are calculated and used for event selection (jAharonian et al.l 119971 : iKrawczvnski et al 



2006). The event selection cuts are optimized a priori using data taken on the Crab Nebula, scaling 
the background and excess rates to account for a weaker source. These selection criteria are termed 
the "standard cuts" . Since 3C 66A is possibly very distant and the observed spectrum is expected 
to be soft, a modified "soft cuts" applies a further a priori optimization of increasing the 9 2 cut 
(the angular distance squared from the position of 3C 66A and the reconstructed shower direction) 
and decreasing the size cut (the number of photo-electrons in an event) for sources with a soft 
spectrum (see Table [I]). Unless stated otherwise, the "soft cuts" were used to generate the results 
presented in this paper. 



The reflected-region model (jBerge et al.l 120071 ) is used for background subtraction. The total 
number of events in the on-source region is then compared to the total number of events in the 
more numerous off-source regions, scaled by the ratio (called a) of the solid angles, to produce a 
final excess. 



4. VERITAS Results 



An excess of 1791 events is observed from the direction of 3C 66A (7257 on events, 31201 off 
events with an off-source normalization ratio a of 0.1752 ) . The excess corresponds to a statistical 
significance of 21.1 a using Equation 17 from iLi Mai ( 19831 ). The distribution of 6 2 is shown 
in Figure [TJ The shape of the excess is consistent with that expected from a point source (68% 
containment in 0.16°). Since lower-energy events are more poorly reconstructed and this data set 
is dominated by low-energy events, the 6 2 distribution is wider than that of a harder-spectrum 
source. Table Q] details the results using soft cuts as well as standard cuts. Note that the use of 
two different sets of cuts has an associated trials factor of two, which has negligible impact on the 
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Table 1: Selection cuts applied to the data for the soft cuts as well as the standard cuts a . Also 
shown are the results of the VERITAS observations of 3C 66A using both the standard analysis 
cuts and soft spectrum cuts. 



Cuts 


e 2 


Size b 




On d 


Off 3 


Alpha f 


Excess g 


Sig. 




[deg 2 ] 


[DC] 


[GeV] 










M 


Soft 


< 0.020 


> 200 


120 


7257 


31201 


0.1752 


1791 


21.1 


Std. 


< 0.013 


> 400 


170 


1258 


5282 


0.1400 


518 


16.0 



"All other cuts outside of those shown here are the same 

6 The photo-electron to digital count ratio is approximately four 

c Post-cuts energy threshold derived using a spectral index of 4.1 

d Number of on-source events passing cuts 

e Number of off-source events passing cuts 

■'Normalization for the off-source events 

9 Observed Excess 
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Fig. 1. — The distribution of 6 2 for on-source events (points) and normalized off-source events 
(shaded region) from observations of 3C 66A. The dashed line represents the cut on 9 2 applied to 
the data. 
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Fig. 2. — The energy spectrum of 3C 66A shown as solid points. The spectrum is well fitted by a 
power law with index T = 4.1 ± 0A sta t ± 0.6 sys (s olid line) . The shaded are a outlines the systematic 
error in the spectral index. Using the models of iFranceschini et al.1 (|2008l ) and assuming a redshift 
of z = 0.444, the de-absorbed spectral index is calculated to be 1.1 ±0.4 showing that the very steep 
measured spectrum could be due to the distance of 3C 66A. This de-absorbed spe ctrum is shown a s 
a dashed-dotted line and points. The MAGIC spectrum with index T = 3.1 from lAliu et al.l ([2009) 
is shown as a dotted line. The Crab Nebula's spectrum divided by 10 is also shown for comparison 
(dashed line). 
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Fig. 3. — Light curve binned by dark period (the time between two full moons) for the full data 
set from 3C 66A. These data indicate night-by-night variability for the second dark period (MJD 
54734 through 54749) but not within any of the other dark periods. The inset details this dark 
period binned by night when the flare occurred. 
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Fig. 4. — Smoothed significance map of 3C 66A. The location of 3C 66A is shown as an open star 
and 3C 66B as the closed star. The cross is the fit to the excess VHE emission resulting in a 
localization of 2h 22m 41.6s ± 1.7s ± 6.0s, 43° 02' 35.5" ± 21" ± 1'30". These data strongly favor 
3C 66A as the source of the 7-ray emission at a significance level of 4.3 a. Note that the bins are 
highly correlated due to an integration over angular space. 
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significance. 

Using the soft cuts, the differential energy spectrum over the energy range ~ 200 GeV to ~ 500 
GeV is determined and is shown in Figure 2. The best fit of a power-law to these data yields an 
index T = 4.1 ± 0.4 

s tat i 0.6 S y S with a chi-square of 1.94 for 2 degrees of freedom. An alternative 
analysis chain confirms this very soft spectrum. Note that there were 1431 excess events detected 
during the flaring period from MJD 54740 through MJD 54749 which accounts for 80% of the 
total. Thus, while the spectrum calculated here is for the full data set, it is dominated by the flare. 
Assuming this power-law spectrum, the observed integral flux for the full data set above 200 GeV is 
(1.3 ±0.1) x lO" 11 cm" 2 s" 1 (6% of the Crab Nebula's flux). By comparison, the 2007-2008 season 
yielded a significance of 2.6 a at a lower average flux above 200 GeV of (3.9 ± 1.6) x 10 -12 photons 
cm -2 s j which is 26% of the flux seen in 2008-2009. Figure [3] shows the integral flux above 200 
GeV from 3C 66A for each dark period (the time between two full moons). The highest flux seen 
from 3C 66A occurred on MJD 54749. Significant variability is seen only during the dark period 
spanning September 25 through October 10 (shown in the inset in Figure [3]), with a chi-squared 
probability of 0.009% for a fit to a constant flux. No statistically significant evidence for variability 
is seen within any of the individual nights. Fits of a constant to the nightly flux in any other dark 
period do not yield a chi-squared probability less than 10%. 

The radio galaxy 3C 66B lies in the s ame field of view as 3C 66A at a separation of 0.12° and is 



also a plausible source of VHE radiation (jTavecchio &: Gl 



isellini 



2008). With the recent detection 



of VHE emission from the 3C66A/B region by MAGIC (jAliu et al.l 120091 ) favoring 3C66B as the 
source and excluding 3C 66A at an 85% confidence level, it is important to determine which of 
these objects is the source of the emission reported here. Thus, a 2-dimensional Gaussian shape 
was fit to the uncorrelated excess of 7 rays, yielding a position of 2h 22m 41.6s ± 1.7s, 43° 02' 
35.5" ± 21", with an additional systematic angular uncertainty of 90". The systematic error has 
been confirmed via optical pointing monitors which are mounted to each telescope. This rules out 
3C 66B as the source of the VHE emission reported here at a significance level of 4.3 a. 3C 66A 
lies 0.01° from the fit position while 3C 66B lies 0.13° away and the total error on the fit is 0.03° 
(see Figured]). In addition to fitting the full data set, fits were made to the data divided into 
high (> 300 GeV) and low (< 300 GeV) energy bands under the assumption that the high-energy 
emission might originate fro m 3C 66B while the low energy emission might come from 3C 66A (see 



Tavecchio Ghisellinil 120081 ) . The fit to the position did not deviate from the measurement using 
the full data set. Correlated variability studies utilizing optical, VHE and HE (30 MeV - 100 GeV) 
bands are underway to verify 3C 66A as the source of VHE 7 rays and will be the subject of a future 
paper. Further restricting the VERITAS data to observations contemporaneous with MAGIC in 
September to December 2007, we calculate an upper limit, assuming the reported MAGIC spectrum 
of T = 3.10, on the flux above 300 GeV from 3C66B to be 1.8 x 10~ 12 photons cm" 2 s" 1 at the 
99% confidence level based on ~ 5 hours of data. MAGIC reported an integral flux based on a 
~ 50 hour exposure above 150 GeV of (7.3 ± 1.5) x 10~ 12 photons cm -2 s _1 for their full data set 
(approximately 1.7 x 10~ 12 photons cm -2 s _1 above 300 GeV). Unfortunately, it is not possible to 
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calculate a spectrum from the 2007-2008 season data due to low statistics. Although VERITAS 
is more sensitive (~ 2x) than MAGIC the brief (~ 5 hour) exposure on the 3C66A/B region in 
2007 does not enable a clear determination of which object was the source of VHE emission in 
2007. However, based on the MAGIC flux, VERITAS expects ~ 700 excess events at the location 
of 3C 66B in the full data set, whereas only ~ 300 are detected and the latter is consistent with 
expectations for spill over from 3C 66A due to the VERITAS point spread function. Therefore, 
if the MAGIC claims of VHE emission from 3C 66B are correct, it must have been considerably 
brighter in 2007 than 2008, and similarly 3C 66A must have been considerably brighter in 2008 
than 2007. 



5. Summary &: Conclusion 



VERITAS has observed the IBL 3C 66A for a total of 32.8 hours good-quality live time from 
September 2007 through November 2008, resulting in the detection of VHE 7 rays with a statistical 
significance of 21.1 a. The average integral flux above 200 GeV is (1.3 ±0.1) x 10" 11 cm" 2 s _1 (6% 
of the Crab Nebula's flux). The differential energy spectrum is well fit by a soft power law with 
index r = 4.1 ± QA stat ± 0.6 sys between 200 and 500 GeV. 

It is thought that the redshift of 3C 66A is z = 0.444 but this measurement is based upon a 
single, poorly detected line ([Miller et al.lll978l ). A definitive measurement of the redshift is needed 
to determine the intrinsic spectrum of 3C 66A, corrected for EBL absorption. The extreme distance 
of 3C 66A, if true, will allow modelers to probe the evolution of the EBL with redshift. Assuming 
the current redshift measurement of z = 0.444, we calculate a de-absorbed spectrum based on 



the EBL models of iFranceschini et al.l (|2008l ) which are based upon recent measurements from the 
optical to the sub- millimeter. The original spectrum along with the corrected spectrum can be 
seen in Figure El While this is not a definitive calculation of the intrinsic spectrum, due to the 
uncertainties in the redshift measurement and in the modeling of the EBL, it illustrates that the 
steepness of the measured spectrum could be due to the distance of 3C 66A. 

The initial announcement of a detection of VHE emission fro m 3C 66A rtswordvhoosh prompte d 
several other groups and instruments to also observe this object ( Tosti 20081 : Larionov et al. 2008 ). 
In addition, the Fermi Gamma-ray Space Telescope detected 3C 66A at a level higher than that 
reported by EGRET. The Swift observatory also monitored 3C 66A over this time period, in the 
X-ray and UV bands. A future paper by the Fermi collaboration, VERITAS collaboration and 
multi-wavelength partners will describe these results and provide details on correlated variability 
as well as a broadband spectral energy distribution for the 2008 data set. 
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